Origins of elastic properties in ordered nanocomposites 
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We predict a diblock copolymer melt in the lamellar phase with added spherical nanoparticles that 
have an affinity for one block to have a lower tensile modulus than a pure diblock copolymer system. 
This weakening is due to the swelling of the lamellar domain by nanoparticles and the displacement 
of polymer by elastically inert fillers. Despite the overall decrease in the tensile modulus of a 
polydomain sample, the shear modulus for a single domain increases dramatically. 
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Polymer nanocomposites are being extensively inves- 
tigated because of the improvement in material prop- 
erties that result from the addition of nanoscopic filler 
particles to the polymer matrix Q, H, H, In addi- 
tion to their practical importance, such composites offer 
diverse scientific challenges, combining ideas from col- 
loid science, polymer physics and chemistry, as well as 
material science. Polymer nanocomposites become even 
more interesting when the polymer matrix consists of a 
block copolymer, capable of self-assembling into a wide 
range of ordered nanoscaled structures — nanoparticles 
can then be sequestered in certain domains to form or- 
dered nanocomposites 0, IE 13 • The simultaneous am- 
phiphilic and colloidal self-assembly taking place in such 
ordered nanocomposites gives them complex structures 
and makes the structure-property relationship par- 
ticularly intriguing. Since there is little understanding of 
the mechanical properties that arise in ordered nanocom- 
posites, we present in this theoretical work a first inves- 
tigation of the origins of the elastic properties of an or- 
dered nanocomposite with spherical nanofillers. 

Buxton and Balazs |0| have studied a phenomenologi- 
cal model of nanosphere filled block copolymer systems in 
which a hybrid Cahn-Hilliard/Brownian dynamics simu- 
lation is used as input to a lattice spring model of the 
elastic moduli. Their approach provides a versatile and 
useful method of predicting properties, but lacks poly- 
meric detail in the elasticity portion of the simulation. 
Furthermore, they examine filled block copolymer sys- 
tems in the solid state, where all morphological evolution 
is disregarded as the system is distorted. 

We examine the elastic properties of a melt state or- 
dered nanocomposite using self-consistent field theory 
(SCFT). SCFT is a coarse-grained, first principles ap- 
proach that been successful in dealing with block copoly- 
mer structure In the framework of this theory, lo- 
cal monomer density profiles of different block copoly- 
mer chemical species are represented self-consistently us- 
ing chemical potential fields. Both the densities and the 
fields are then used to determine the free energy for the 
system, and if desired, the internal energies and entropies 
can be explicitly calculated. SCFT has been extended to 



deal with hard nanosphere/block copolymer nanocom- 
posites by the incorporation of a density functional the- 
ory particle contribution 0, Further, Tyler and 
Morse have demonstrated that the linear elastic behav- 
ior of a melt block copolymer system, which is quasi- 
statically deformed can be well characterized using SCFT 
[l^ . We have recently adapted this approach to an effi- 
cient real space, pseudo-spectral method 15j and found 
an increasing elastic modulus in multiblock copolymer 
systems as a function of block number, in qualitative 
agreement with experiment [l^ . Here, we combine these 
two advances in order not only to predict the effect on 
the elastic properties of adding nanoparticles to a block 
copolymer melt, but also to explain the physical origins 
of the observed effects. 

We will study the prototypical system consisting of a 
symmetric AB diblock copolymer melt in the lamellar 
phase with added spherical nanoparticles that have an 
affinity for the A block of the copolymer. Consequently, 
a lamellar morphology with the particles sequestered in 
the A phase is being considered, and the system's tetrag- 
onal symmetry is elastically characterized by just five 
independent non-zero components of the elastic modu- 
lus tensor. Additionally, the system is in a melt state so 
that deformations parallel to the lamellar structure have 
no effect on the free energy of the system |0. We are 
thus left with only two relevant moduli, K33 and K44. 
Therefore we deform the system quasi-statically in two 
ways; it is subjected to an extension/compression, and 
to a simple shear. These deformations allow us to deter- 
mine the K33 and K44 components of the elastic modulus 
tensor jl7| . corresponding to extension/compression and 
shear moduli, respectively. These components are found 
by taking the second derivative of the SCFT free energy 
with respect to the relative deformation. Greater detail 
on the methodology can be found in Ref. [isj . 

We chose a system with a segregation of = 25 be- 
tween the A and B blocks, with the particles considered 
to be of the A species, x is the Flory-Huggins monomer 
segregation parameter and N is the degree of polymeriza- 
tion of the entire diblock. The particles radius was cho- 
sen to be 0.725i?g, where Rg is the unperturbed radius 
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TABLE I: Components of the _R'33 elastic modulus for 0% 
and 15% added nanoparticles. The change in the modulus 
between the pure diblock and the filled diblock system is 
recorded in the last column. 
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TABLE IL Components of the K44 elastic modulus for 0% and 
15% added nanoparticles. The change in modulus bet-ween 
the pure diblock and the filled diblock system is recorded in 
the last column. 



of gyration of a diblock molecule. Finally, the particle- 
to-diblock volume ratio was ~ 3.6, and a 15% volume 
fraction of spherical fillers was added. The system was 
deformed in the two ways described above and compared 
with a neat diblock system similarly deformed. The K^^ 
and K44 moduli in each case were used to find a tensile 
modulus by averaging over a polydomain sample accord- 
ing to the Hill prescription as described in Ref. 0| . 

For the pure diblock system a dimensionless tensile 
modulus of ^ 0.4 was found while for the nanocompos- 
ite system, the modulus was ~ 0.3. Thus we predict 
that the addition of nanoparticles will weaken the ma- 
terial. Normally, one expects the addition of nanofillers 
to strengthen the composite, rather than weaken it 0, 
although a reduced modulus has been observed experi- 
mentallyin exfoliated layered silicate/triblock nanocom- 
posites |3- The reduced modulus found for our present 
system can be explained by decomposing the K33 and 
K44 moduli as shown in tables ^ and These tables 
show total K33 and K44 moduli values, respectively, for 
0% and 15% added fillers, as well as the difference be- 
tween the filled and unfilled system moduli. The tables 
also give the contributions to the moduli of the internal 
energy (C/), translational entropy {St), and the A and 
B conformational entropies of the diblocks (5*^ and Sb)- 
Also listed are the ideal gas (Sid) and steric {Sst) con- 
tributions of the filler particles. For information on how 



these components are calculated, see Ref. 0. Compar- 
ing the total and K44, it is observed that the main 
contribution to the tensile modulus arises from the ex- 
tensional modulus Furthermore, the main drop in 

K33 from the 0% to the 15% case is a result of the inter- 
nal energy contribution and the A block conformational 
entropy contribution. 

These two contributions can be examined separately. 
Upon extension/compression, SOFT shows that the 
lamellar interfacial width remains practically unchanged 
in both the filled and unfilled cases. Furthermore, the A 
and B domains are well segregated both before and after 
deformation, indicating that there are few A monomers in 
the B region and vice versa. Thus the absolute amount of 
energetically unfavorable AB contacts is the same before 
and after deformation, whereas the domain size changes. 
Adapting an expression of Matsen and Bates y^j, the 
internal energy contribution to the free energy is 
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where T is the temperature and ks is Boltzmann's con- 
stant. (pa{^), <PB(r) and (pp{r) are the local volume frac- 
tions of A and B monomers, and particles, respectively. 
The integral in Eq. JQ) is constant under the conditions 
described above, so that the internal energy is inversely 
proportional to the volume. The lamellar morphology is 
one dimensional, so that the internal energy can be writ- 
ten as inversely proportional to the equilibrium domain 
size d*. 
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where a is a constant. A relative deformation e that 
leaves the interfacial width and the bulk mixing un- 
changed, changes Eq. ^ into 
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The SCFT internal energy for extensions/compressions is 
shown in Fig. QJa). Both for pure diblock and for filled 
systems, the behavior reflected in Eq. ^ is observed. An 
estimate of the modulus K^^ can be found by taking the 
second derivative of Eq. Q with respect to the relative 
deformation, which gives 
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showing that the U contribution to the if 33 modulus is 
also inversely proportional to the equilibrium spacing. 
SCFT shows that the addition of filler particles enlarges 
the equilibrium domain size of the nanocomposite com- 
pared to the pure diblock system. Equation Q then 
indicates that the modulus will drop, as observed. In 
other words, the modulus is weakened partially because 
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FIG. 1: Internal energy contribution to the free energy versus 
relative distortion for extension/compression (a) and shear 
(b) of a diblock system. The neat diblock free energy is plot- 
ted with the solid line whereas the 15% filled system is shown 
with dashed line. In (a), negative deformations are compres- 
sions while positive deformations are extensions. In all plots, 
energies have been zeroed around the equilibrium spacing, 
which is represented by e = 0. 
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FIG. 2: Conformational entropy contribution to the free en- 
ergy versus distortion for extension/compression (a) and shear 
(b) of a diblock system. The neat diblock free energy is plot- 
ted with the solid lines whereas the 15% filled system is shown 
with dashed lines. In (a), negative deformations are com- 
pressions while positive deformations are extensions; the Sa 
contribution is shown in (a). The Sb contribution is shown 
in (b). In all plots, energies have been zeroed around the 
equilibrium spacing, which is represented by e = 0. 



there is less interface per volume in the nanocomposite 
compared to the pure diblock. 

The A block configurational entropy also contributes 
to the overall drop in modulus. The filler particles have 
no configurational entropy; in the filled system the Sa 
energy portion will rise (drop) under extension (compres- 
sion) at a slower rate because there is a smaller fraction 
of chains to stretch (relax). This can be seen in Fig. Ufa). 
Thus filler particles weaken the material because they dis- 
place polymers that have stretching energy, which could 
contribute positively to the elastic modulus. In other 
words, the modulus is weakened partially because there 
is less polymer and more elastically inert filler per vol- 
ume. 

As mentioned, the shear modulus K^^ only makes a 
small contribution to the tensile modulus, but it is in- 
teresting to examine it nonetheless. Table ^ shows that 
while the if 33 modulus drops upon the addition of fillers, 
K44 increases by more than 400%. From tabled this 
increase can be seen to be due to a very large increase 
in the B block conformational entropy contribution to 
the modulus. This is partially reduced by a drop in the 
internal energy contribution, but is still very large. 

As in the extension/compression case, the interfacial 
profiles do not change significantly upon shearing of the 
sample. The consequence is that the internal energy 
again is inversely proportional to the volume of the sys- 
tem, as previously explained. The B block conforma- 



tional entropy contribution to the free energy behaves 
similarly; from Matsen and Bates the B block en- 
ergy is 
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where pj{r) is the distribution of diblock junction points, 
wb(j) is the chemical potential field for the B monomer 
distribution, and q^{r,s) is a SCFT propagator. A de- 
tailed explanation of Eq. (jSJ can be found elsewhere 
[nl nil . If we for the moment ignore the contribution of 
the integral in Eq. (|SJ), the B conformational contribu- 
tion can be written as 
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where /3 is a constant. A sin 6' has been added to the 
denominator of Eq. © to account for the effect of a 
shear change in volume. When the sample is sheared, the 
volume is reduced; is a measure of the amount of shear 
with 9 = Ti 12 representing a non-sheared system. See 
Ref. 'l5"| for more explanations. The internal energy (0) 
will be similarly affected. From Ref. Ill, 

the shear angle 

is related to the relative distortion through e = cot 6. The 
internal energy and B conformational entropy are then 
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Figures n^b) andj^b) show that the internal energy and 
B block contribution to the free energy obey the rela- 
tionships and JHJ, respectively, for both the filled and 
unfilled systems. The moduli for these are found through 
the second derivative with respect to e and are 
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Both K^^ and K^^ are inversely proportional to the equi- 
librium volume, with the important difference of a minus 
sign. Thus as the equilibrium volume is enlarged upon 
the addition of fillers, the Sb contribution to the modu- 
lus becomes a smaller negative number, so the material 
is stronger. This effect is somewhat reduced as the U 
contribution to the modulus will become smaller for the 
larger equilibrium volume of the filled system. In other 
words, diblock molecules' entropies want to help shear 
the system (see Fig. EJb)). The larger domain size of 
the filled system means that there is less stretching en- 
ergy per volume (see eq. ((HI), the molecules are more 
relaxed to begin with at equilibrium, and so are less in- 
clined to help deform the system, which makes the ma- 
terial stronger. This effect is reduced since there is less 
interface per volume in the filled system. The A block 
is not treated the same way, since with the fillers added 
the Sa contribution to the free energy cannot be written 
in a form such as Eq.l^. 

In summary, we have calculated the tensile modulus 
for a neat diblock copolymer system and for a diblock 
nanocomposite with 15% added nanospheres. Both sys- 
tems were considered to be in the lamellar phase, and 
deformations were applied quasi-statically. The elastic 
modulus tensor components and K44 were found and 
used to derive a tensile modulus for a polydomain sample. 
It was found that the addition of nanoparticles weakened 
the material. This was attributed to the larger lamellar 
domain size of the equilibrium filled system — it had less 
interface per volume — and to the displacement of poly- 
mer by the filler particles — there was less elastic polymer 
per volume. Although, the shear modulus contributed 
little, it was acted on by similar mechanisms, with the 
result that it increased dramatically upon the addition 
of nanofiUers. Given that ordered nanocomposites of the 
sort described here can now be realized, we believe our 
predictions should be amenable to experimental verifica- 
tion. It would be interesting to change the distribution 
of the particles in the block copolymer through changes 
of wetting properties, particle size, or volume fraction. 



Such changes in structure could significantly change the 
properties. 
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